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Abstract 

We compute the static potential of adjoint sources in SU(2) Yang-Mills theory in 
four dimensions by numerical Monte Carlo simulations. Following a recent calcu- 
lation in 2+1 dimensions, we employ a variational approach involving string and 
gluelump operators and obtain clear evidence for string breaking and the satu- 
ration of the potential at large distances. For the string breaking scale we find 
~ 1.25fm, 2.3ro, or in units of the lightest glueball, r(,mo++ ~ 9.7. We fur- 
thermore resolve the first excitation of the flux-tube and observe its breaking as 
well. The result for r?, is in remarkable quantitative agreement with the three- 
dimensional one. 
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1 Introduction 



The most essential property of QCD is confinement: the energy of a static quark- 
antiquark QQ pair grows linearly with their separation r, because the colour flux between 
the two charges is negligible outside a thin tube, the string, whose energy grows in pro- 
portion to its length. This remains true up to a distance where the string breaks: the 
energy becomes sufficient to create a pair of mesons, Qq and Qq, which can be taken 
arbitrarily far apart at no extra energy cost. The static potential thus levels off at a 
value 2E(Qq). Similarly, the force between a pair of adjoint sources Q a dj is screened by 
gluons g, so that the adjoint string breaks, even in the absence of quarks, when a pair of 
"gluelumps" Qadjd can be created. In numerical simulations of lattice QCD it has been 
far easier to observe the string itself than its breaking. In spite of extensive searches for 
colour screening |IJ, the evidence is still preliminary at best 0. 

The essential reason lies in the projection properties of the Wilson loop, whose change 
from an area law to a perimeter law behaviour would signal string breaking. The difficulty 
with seeing the perimeter law in simulations stems from the fact that at r > r& the state 
with unbroken string between the sources is still an eigenstate of the Hamiltonian, albeit 
an excited one. By construction, the Wilson loop has good overlap with this state, which 
remains the case also for r > r&. If it also has poor overlap with the two-gluelump 
state, then for r > r& the latter, which is the true ground state, will only dominate the 
correlation function for very large correlation times t, which according to the definition of 
the potential should be taken to infinity for the loop to yield the ground state. In most 
practical simulations, the limit t — > oo is not realized, but t is typically a few lattice 
spacings when the signal is lost in noise. If this happens before t is large enough for the 
ground state to dominate the correlation function, the ground state will be missed and 
an excited state will be extracted instead, leading to the apparent conclusion that the 
string does not break ||. 

A way to circumvent this difficulty is to supplement the calculation with an addi- 
tional operator having good projection onto the two-gluelump state. This enables an 
unambiguous identification of the ground state also for r > r&, even if the perimeter 
behaviour of the Wilson loop is hidden in noise. After a first attempt for the adjoint 
potential || , this procedure has been first sucessful for the fundamental static potential 
of the SU(2) Higgs model in its confining phase in (2+1) as well as in (3+1) || di- 
mensions, where it clearly demonstrated the flattening of the potential as well as mixing 
between the Wilson loop and the broken string state in a very narrow region around 
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r&. The same findings were reported from simulations of the adjoint potential in (2+1) 
dimensional Yang- Mills theory [[?]. Q . 

In this work, we study string breaking in the adjoint potential of SU(2) Yang-Mills 
theory in four dimensions. In practical terms, the advantage over QCD is that no dy- 
namical quarks need be considered. The price to pay is that the adjoint string tension 
is larger than the fundamental one, approximately by the ratio 8/3 of the Casimirs, so 
that Wilson loops obeying an area law disappear more quickly in the statistical noise, 
which turned out to hamper the earlier study {§. 

Section 2 describes our calculational method for on- and off-axis correlations. Section 
3 presents our simulations and Section 4 our numerical results. Section 5 contains a 
summary and outlook. 



2 Calculational method 



We work with the Wilson action 

S[U] = 



X )1<V 



(1) 



where denotes a plaquette of links in the fundamental representation, f3c = 4/g 2 
and g is the bare gauge coupling. The links in the adjoint representation, Af^, are related 
to the fundamental U^s by 



where a a are the Pauli matrices. 



^Tr (a a U,(x)a b Ul(x 



(2) 



The standard operators to extract the static potential are rectangular Wilson loops 
with side lengths r and t in a spatial and the time direction, respectively. Physically, a 
Wilson line S(x, y) is interpreted as colour flux propagating along that line from x to y. 
Accordingly, a straight Wilson line in t-direction corresponds to a static colour source, 
whereas a Wilson line in r-direction may be used to describe a flux tube between x and y. 
A rectangular Wilson loop then represents the correlation function of a string of length 
r over a time interval t. To consider the case of adjoint sources, an adjoint Wilson loop 
is needed, which takes the form 



where 



G S s(r,t) = W^(r,t)=\W(r,t) 



W{r, t) = Tr \S{0, rj)%j, rj + t6)S\t6, rj + t6)S f (0, t6) 



(3) 
(4) 
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is the fundamental Wilson loop, S(x, y) is a straight line of fundamental links connecting 
the sites x and y, and j and denote unit vectors in a spatial and the time direction 
respectively. The static potential in the adjoint representation is defined in terms of the 
exponential decay of the adjoint Wilson loop, 



V(r) 
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Km -ln[W adj (r,t)]. 



(5) 



In the region of linear confinement the Wilson loop obeys an area law, whereas for 
distances beyond some breaking scale a perimeter law is expected. This perimeter 
behaviour for large separations r, however, has not been observed to date. 

The correlation function for a bound state of a static adjoint colour source and a 
gluon field, called gluelump in the literature || , is given by the non-local gauge-invariant 
operator 



G G {t) 



Ti(P(x)a a )SZ b dj (x, x + t6)Tr(P t (a; + tb)a h 



Tr 



P(x)S(x, x + t6) (P\x + t6) - P(x + tO)) S^(x, x + tO) 



with the adjoint representation Wilson line 



S%(x,y) = -Tr (a a S(x,y)<r b S\x,yj) . 



(6) 



(7) 



Here P(x) denotes the "clover- leaf" of all four plaquettes with the same orientation which 
emanate from the endpoint x of the adjoint Wilson line, summed over the three spatial 
planes || . 

We now follow the procedure proposed in M and applied in three dimensions 
construct an operator projecting on two of these bound states at distance r by 



to 



G GG (r,t) 



Tr 

x Tr 



p(o)s(o,to) (p ] (to) - P(to)) 5 f (0,t6) 

P(rj)S(rj, r] + tO) (P f (rj + tO) - P(rj + tO)) S\rj, rj + tO) 



Finally, correlations between a string and a gluelump state, and vice versa, may be 
described by 



G SG (r,t) 
G GS (r,t) 



Tr 



Pt(t6) - P(tO)) S\0,tO) S(0,rj)S(r),r) + tO) 



x 



P(rj + t6)^ f (rj, rj + t6)S t (0, rj)S(0, tO) 



Tr 



P f (0) - P(0)) S(0, tO) S(t6, rj + t6)S\rl rj + tO) 



x 



P(rj)S(rj, rj + t6)S\tO, rj + tQ)S ] (0, tO) 



(9) 
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The static potential, its excitations and the mixing between gauge string and two- 
gluelump state can then be extracted from measurements of the matrix correlator 

V G Gs{r,t) G GG {r,t) J 

We also keep the single gluelump operator G G from in our simulations, in order to 
check whether G GG defined in (||) indeed has a good projection onto a two-gluelump 
state, for which one expects E GG ~ 2E G . 

A particular difficulty in observing string breaking comes from its suddenness. As 
the distance between the sources increases, the ground state switches abruptly from an 
unbroken string overlapping almost exclusively with the Wilson loop, to a broken string 
overlapping almost exclusively with two gluelumps. Mixing among the two states, which 
proves that one passes continuously from one ground state to the other, occurs only if the 
two states are very close in mass. This corresponds to a very narrow window of distances, 
which becomes elusive or entirely unobservable on a coarse lattice M . For this reason, we 
supplement our study of on-axis correlations as above with that of off-axis correlations, 
which offer far more discrete Euclidean distances. To avoid computing spatial parallel 
transporters from x to y as in eqs. (§,|8],|9|) for off-axis separations, we instead fix the 
gauge and measure in a fixed gauge the generic adjoint correlator between two time-like 
"segments" 

G fix (r,t) = (| TrS(0,t6)S t (rj,rj + t6) | 2 -1). (11) 

Any gauge which preserves the transfer matrix is acceptable. The choice of gauge will 
determine the projection of the above correlator onto the string or the two-gluelump 
state. To project onto the broken string, Gfi X ~ G GG , we simply choose the gauge 
(defined up to an irrelevant Abelian rotation exp 

(i0cr 3 )) which brings F = ±(P - pt) to 
the colour direction <7 3 at every point: F 1 (x) = F 2 (x) = Vx. The required local gauge 
transformation can be determined and applied independently at each lattice site. To 
project onto the unbroken string, Gfi X ~ Gss, we need a smooth gauge and choose the 
adjoint Coulomb gauge. Correspondingly, for the mixed elements G$ G the gauges on the 
two timeslices in the correlator G^ x are different. In order to avoid possible systematic 
effects from lattice "Gribov" copies, we discard the usual implementation of Coulomb 
gauge on the lattice by iterative local extremization |1(J in favour of the unambiguous 



prescription of Vink and Wiese |TT[. Gauge-fixing is accomplished by rotating the two 
lowest-lying eigenvectors of the adjoint covariant 3d Laplacian to direction cr 3 and into 
the plane (01,03) respectively in colour space. This gauge condition corresponds to the 



Laplacian Center Gauge described in ||12|| , with only the three spatial directions used for 
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the Laplacian. It reduces to Coulomb gauge for the adjoint field in the continuum limit. 
It is unambiguous except when either of the lowest two eigenvalues of the Laplacian 
becomes degenerate, which signals a genuine (not a lattice artifact) Gribov copy. Such 
a degeneracy did not occur over the whole course of our simulations. The computation 
of the lowest-lying Laplacian eigenvectors is performed with the use of the ARPACK 
package JT3|] . Note that it requires no more computing effort than the usual iterative 
procedure. 

Nevertheless, gauge-fixing remains expensive. In principle, we would like to perform a 
variational study over a variety of gauges. To keep computer resources within our budget, 
we had to limit ourselves to the two gauges above, where the links used to construct F 
or the covariant Laplacian have been smeared 10, respectively 20 times. The price we 
pay for such a truncation of the variational basis is a poorer determination of the ground 
state. Using smeared rectangular Wilson loops yields more accurate results at cheaper 
cost. The specific purpose here is to exhibit the progressive mixing of the two variational 
states with distance. 

The energy levels extracted from the matrix correlator G(r, t) in eqn. (|T0|) are linearly 
divergent with decreasing lattice spacing and hence do not have a continuum limit. This 
divergence is due to the self-energy of the static sources which, although perturbatively 
computable, cannot be absorbed by renormalization into a parameter of the theory. It 
reflects the fact that the static potential is defined only up to an arbitrary (infinite) 
constant, and does not itself constitute a finite physical quantity. On the other hand, the 
confining force and the string breaking scale are defined by energy differences such that 
the divergence cancels out, and these quantities do have a finite and physical continuum 
limit. 



3 Simulation and analysis 

In order to improve the projection properties of our operators, we employed the standard 



smearing algorithm [TJ[] to obtain smeared spatial link variables of unit length, which 
were then used instead of the original ones in constructing the correlation functions 
defined above. All links in the time directions were left unsmeared such that the transfer 
matrix remains unaffected by our smearing procedure. As a basis of operators used in 
the matrix correlator eq. ([TOD, we chose two different link smearing levels for the spatial 
Wilson lines and one smearing level with good projection for the clover-leaves P. Thus, 
our correlator G(r, t) represents a 3 x 3 matrix for the gauge-invariant calculation. In 
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the case of the off-axis, gauge-fixed correlators we kept only one smearing level for the 
unbroken string, giving a 2 x 2 matrix for G(r,t). 

The procedure we follow to diagonalize G(r,t) by means of a variational calculation 
has been described in detail in the literature |T3|, [Tj| , and its application to the calculation 
of the adjoint potential has already been discussed || ^, [7], ||. Here we use the same 
algorithm as in M and for completeness briefly repeat the outline of the procedure. 
Each element Gy of the matrix correlator represents a correlation function which can be 
written as 

G«(r,t) = (0 i (t)0 i (O)> l (12) 

where the 0j's represent a spatial gauge string of length r at a given smearing level 
for i = 1,2, and a two-gluelump operator at a different smearing level for i = 3. In a 
temporal gauge, where all time-like links are unity, this expression describes indeed the 
measured correlation functions. Since all Gij are manifestly gauge invariant, this remains 
true without gauge fixing. The variational diagonalization of G(r, t) consists of solving 
the generalized eigenvalue problem || [T7| 

G(r, t) v% t ) = \(t, t )G(r, *„)«*(<, t ), t> t . (13) 

From the eigenvectors v 1 one may construct the corresponding eigenstates 

$i = Ci 53 V l k (fik = 53 a ik4>k, (14) 
k k 

which are superpositions of the string and gluelump operators used in the simulation. The 
constants q are chosen such that $j is normalized to unity. The diagonalized correlation 
matrix may then be written as 

3 3 

Ti(r,t) = ($i(t)$i(0)) = J] aijUih 4)^)^(0) = 53 aijdikGjkit), (15) 

j,k=i j,k = i 

and represents the (approximate) correlation functions of the eigenstates of the Hamil- 

tonian. To check the stability of the procedure we have performed the same calculation 

for to = 0, t = a, 2a and obtained consistent results in both cases. The coefficients 

take values between zero and one and quantify the overlap of each individual correlator 

Gik with the correlators of the mass eigenstates, IV They characterise the composition of 

the eigenstates in terms of the original operators as well as the degree of mixing between 

operators. 

We extract the effective masses of the eigenstates $j by assuming a single exponential 
decay of their correlation functions, 

(16) 



K ff {t) = - In 



Tj(r,t+1) 

r«(r,t) 
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The effective masses in lattice units are so heavy that the observed "plateaux" only 
extend over two or three t- values before the signal is lost in noise. In general we therefore 
approximate the z-th eigenstate of the static potential in lattice units by aVi(r/a) = 
m l e ff{2). This leaves room for a systematic error because the asymptotic plateau values 
might be attained only at larger t and thus the energies be overestimated. In the case 
of the two-gluelump operator this error can be ruled out because we find it to be fully 
compatible with twice the single gluelump mass, for which we do see a good plateau. 
Only for the Wilson loop does this source of uncertainty remain. Since string breaking 
occurs when the energies of the broken and unbroken string are equal, overestimating 
the latter may lead to underestimating the string breaking distance r&. 

The size of the operator basis used here is significantly smaller than in the corre- 
sponding three-dimensional study, where up to eight operators were used |§. The reason 
simply is that our present computer resouces do not allow for a larger basis in a four- 
dimensional calculation. In consequence, whereas we may be fairly confident about our 
mass estimates for the ground state potential, the first and especially second excited 
states are likely to be contaminated by higher excitations. In particular, since we have 
only one two-gluelump operator in our basis, we are unable to compute the first excita- 
tion of the gluelump and thus cannot determine the correct ordering of the excited states 
for distances beyond the breaking scale, r 3> r&. We shall come back to this point in the 
discussion of the results. 

In order to examine the scaling behaviour of the potential, we have worked at two 
values of the bare coupling, (5 = 2.4 on a 32 4 lattice, and (3 = 2.5 on a 40 4 lattice. We did 
not find a value for vq(/3 = 2.4) in the literature, and hence based our scaling analysis 
on the expression 

(3\n2 ~ 
~0.205_ 

which describes the scaling of the fundamental representation static potential with the 



a((3) = 400 exp 



fm , (17) 



lattice spacing in the range of interest |L8| . For better comparison with other simulations 
we also quote our results at (3 — 2.5 using ro((3 = 2.5) = 6.39(09)a |19| . 

In three dimensions an explicit test for finite volume effects has been performed, with 
the finding that a spatial size L « 2.4r& is large enough for to be free from finite size 
effects ||. In the present work we shall find that L > 2.8?"b for both lattice spacings and 
we expect no finite size effects. 

At (5 = 2.4 we employed a maximum of 40 smearing steps, whereas at (3 = 2.5, 
60 steps were necessary to reach approximate saturation in the lowering of the effec- 
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tive mass of the Wilson loop. For both lattice spacings, 50 Monte Carlo sweeps were 
inserted between successive measurements, where one heatbath sweep was followed by 
four overrelaxed sweeps. For (3 = 2.4 we collected 400 measurements for the gauge- 
invariant approach and 580 for the gauge-fixed off-axis operators, whereas for (3 = 2.5 
190 gauge-invariant measurements were sufficient to yield comparable statistical errors. 
The latter were estimated using a jackknife procedure, except in the case of the off-axis 
data, where we performed a bootstrap analysis to extract errors more reliably. We have 
checked the decorrelation of our measurements by observing stability of the statistical 
errors when the data were binned. 

4 Numerical results 

As was discussed in Section 2, we measured the mass of the gluelump as a reference to 
check the projection properties of the two-gluelump operator. The correlation function 
G G (t) is relatively easy to measure as it has a good signal and plateaux in the effective 
masses can be identified. For our two lattice spacings we extract 



M G {f3 = 2.4) = 2541(17)MeV, M G {(3 = 2.5) = 2938(23)MeV [r M G = 8.12(13)] . 



In the second line the gluelump mass is given in continuum units. Recall however, that 
there is no physical meaning to these numbers, but the increase on the finer lattice signals 
the linear divergence due to self-interactions of the static source. 

Figure [l] shows the three lowest states of the adjoint static potential in continuum 
units, as obtained from the gauge-invariant calculation at both lattice spacings. We begin 
by discussing the ground state. At small distances we observe the expected linear rise 
corresponding to an area law behaviour of the Wilson loop. The ground state then satu- 
rates rather suddenly at r « 1.25fm, maintaining an energy value of twice the gluelump 
mass, in accord with the expectation of string breaking. This observation is corroborated 
by the analysis of the operator content of the ground state, which is displayed in Figure 
(a). In analogy with previous analyses in (2+1) dimensions || and the fundamental 
potential in Higgs models || ||, we find the Wilson loop to have nearly full projection 
onto the ground state with unbroken string for r < r^, but practically no projection onto 
the saturated ground state for r > r b . 

On the contrary, it maintains its projection onto the unbroken string state which is 
still in the spectrum but now as the first excited state. Conversely, the gluelump operator 



aM G (l3 = 2.4) = 1.54(1) 
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Figure 1: The adjoint potential and its first two excitations for (3 = 2.4 (open symbols) 
and 13 = 2.5 (full symbols). The dashed lines mark the error band for 2M@. 

starts out with some projection onto the ground state at small distances, which may be 
due to gluon exchange between the gluelumps, or simply be an artefact of the smearing 
procedure, since smeared gluelump operators overlap for small distances. However, this 
overlap with the ground state rapidly decreases with distance up to the breaking scale, 
where it projects fully onto the ground state, the latter now being a two-gluelump state. 
The crossover region with appreciable mixing between the two operator types is even 
narrower than was reported in the three-dimensional calculation, presumably because 
we are further from the continuum limit than the simulations in || . 

The ordering of states is completely analogous to the three-dimensional case [Q]. In 
particular, at small distances the first excited state corresponds to an excitation of the 
gauge string, as follows from its linear rise as well as from the overlaps displayed in Figure 
[| (b). It is moreover visible that it receives more contributions from the lightly smeared 
Wilson loop than from the highly smeared Wilson loop. Both figures furthermore show 
quite clearly the breaking of the excited string. In fact, the observed mixing between the 
Wilson loops and two-gluelump operator appears rather more pronounced than in the 
very rapid crossover of the ground state string breaking. We ascribe this to the fact that 
it happens at smaller r, where the Wilson loops are not as large and thus give a better 
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(a): ground state 



(D): tirst excited state 



(c): second excited state 
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Figure 2: The overlaps of the operators in the simulation with the lowest states of 
the potential, as obtained at (5 = 2.5. Open triangles and squares: Gss after 40 and 60 
smearings, respectively; Filled circles: Gqg- 

signal for the variational calculation. After the breaking of the excited string, the first 
excited state of the potential is the two-gluelump state receiving its main contribution 
from the corresponding operator, until this level is attained by the ground state string. 

For r > rfe we do not see saturation of the linearly rising first excited state, in contrast 
to the three-dimensional findings |§ . This is a consequence of our limited basis consisting 
of just one gluelump-pair operator, which prohibits a calculation of excited gluelump 
states. For the same reason, we expect the excited states to be still contaminated by 
higher excitations. 

We now follow the prescription given in || to determine the value of the breaking 
scale r b in continuum units. A suitable definition of r b for this purpose is by the distance 
where the string and the two-gluelump state have equal energy, 



A = Ec - E, 



GG 



0. 



r=r b 



(19) 



To solve this equation we interpolate the energy difference between the string and the two- 
gluelump states to the point where it vanishes. Systematic uncertainties are estimated by 
varying the number of points used for the interpolation. From this procedure we obtain 
the results 



r b /a{(3 = 2.4) = 10.6(8), r b /a{f3 = 2.5) = 14.4(6) 

r b (P = 2.4) = 1.27(8)fm, r 6 (/3 = 2.5) = 1.23(5)fm [r b /r = 2.25(10)] 

which exhibit satisfactory scaling behaviour. 



(20) 



In order to have a higher resolution in the range of string breaking and to obtain a 
better signal for mixing of the operators in the ground state we have performed a gauge- 
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fixed calculation including off-axis correlations at irrational distances in lattice units. 
We have checked explicitly that the on-axis results obtained in this manner agree with 
those from a gauge- invariant calculation with a 2 x 2 matrix. It turns out, however, that 
Wilson loops yield a better signal than the spatial correlation of gauge-fixed sources, 
which need more time-like distance t to attain their asymptotic behaviour. The overlaps 
of the two operators with the ground state in the string breaking region are shown for 
(3 = 2.4 in Figure |3|. Since this is the coarser lattice, the crossover in the overlaps is very 
rapid, even with the higher resolution. Despite their poorer quality, the data show clear 
evidence for progressive mixing between the two operator types, which is absent in the 
gauge-invariant calculation at the same f3, where the ground state correlation is always 
composed entirely of either Wilson loop or two-gluelump operator. Finally, we remark 
that the one would extract from this figure is lower than the one extracted from the 
gauge- invariant calculation in eqn.(^), although they are statistically compatible. We 
nevertheless interpret this difference as a systematic effect rather than a statistical one, 
ascribing it to the severe limitation that a 2 x 2 basis poses. 




9 9.5 10 10.5 11 
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Figure 3: The overlaps of the operators with the ground state of the potential, as 
obtained at (3 — 2.4 in the gauge-fixed calculation. Open squares: Gss! Filled circles: 
Ggg- 
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5 Conclusions and outlook 



We have measured the adjoint static potential in 4- dimensional SU(2) Yang-Mills theory. 
A scaling study at j3 = 2.4 and 2.5 indicates that string breaking occurs at a distance 
rh pa 1.25fm, 2.3r . By measuring off-axis correlations in a fixed gauge, we have observed 
mixing between broken and unbroken string operators, over a distance smaller than one 
lattice spacing (~ 0.1 fm). We have also resolved the breaking of the excited string, 
which occurs at distance ~ 0.75fm, 1.4ro- 

Our findings are in line with the expectation that the string should break where 
the linearly rising part of the potential extracted from pure Wilson loop calculations 
intersects with the two-gluelump state |4| . The suddenness of the breaking, characterized 
by the near-absence of mixing noted above, as well as the ordering of the excitations 
are similar to what happens in (2 + 1) dimensions ||. In fact, the similarity is even 
quantitative: the string-breaking distance in units of the lightest glueball mass p0[ is 



practically the same (r;,m ++ ~ 9.7 versus 10.3(±1.5) in Q). This suggests that identical 
effective mechanisms are at work. 

We find no reason to expect surprises in the case of 577(3). At (3 — 6.0, the gluelump 
mass is known (aMc = 1.33 [^1|]), and string breaking should occur at ~ 12.5a [ pl|j , 
which corresponds again to a similar value of rf,m ++. To observe adjoint string breaking 
at (3 = 6.0, a 32 4 lattice should be marginally sufficient. Such a project is well within 
reach of better-equipped lattice groups. 

The only technical difficulty here is to ascertain that broken and unbroken string 
operators are not completely decoupled, but do mix in the close vicinity of r^. The mixing 
"window" is remarkably narrow, and it would be interesting to find the dynamical reasons 



for this. According to a recent suggestion [f22jl , the suddenness of string-breaking might 
be caused by a qualitative change in the world-sheet spanned by the 2 static charges, 
which tears abruptly as the two sides are taken apart. If this change is of topological 
nature, one may wonder whether mixing should occur at all, or whether the string simply 
"snaps" . 
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